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The stability properties and perturbation-induced dy- 
namics of the full set of stationary states of the nonlin- 
ear Schrodinger equation are investigated numerically in two 
physical contexts: periodic solutions on a ring and confine- 
ment by a harmonic potential. Our comprehensive studies 
emphasize physical interpretations useful to experimentalists. 
Perturbation by stochastic white noise, phase engineering, 
and higher order nonlinearity are considered. We treat both 
attractive and repulsive nonlinearity and illustrate the soliton- 
train nature of the stationary states. 



I. INTRODUCTION 

The one-dimensional nonlinear Schrodinger equation 
(NLSE) is ubiquitious. Among other natural phenom- 
ena, it models the dilute-gas Bose-Einstein condensate 
(BEC) in the quasi-one-dimensional regimdl light pulses 
in optical fibersEl, Bose-mndensed photonsa, helical ex- 
citations pof a vortex lineta, and spin waves in magnetic 
materialsO. In this article application of the NLSE to the 
BEC is emphasized. As currenttjO and proposed! BECs 
are studied in traps of differing topologies, issues of sta- 
bility for periodic solutions on a ring and for confinement 
in a harmonic potential are both-jaf interest. 

In a recent pair of articlesBEJ the full set of peri- 
odic solutions of the stationary NLSE on a ring was 
presented in closed analytic form. Here we extend our 
previous work by numerically studying the stability of 
such stationary states in response to initial and con- 
tinuous stochastic white noise, phase engineering, and 
higher order nonlinearity, and by undertaking a similiar 
set of studies for the harmonic potential. In the BEC 
stochastic white noise arises from the effect of the ther- 
mal cloud of uncondensed atoms on the condensate, and 
is nominally a controllable parameter; simnle phase pro- 
files are already in use in experimentst3t3; and higher 
order nonlinearity models both mea^ti^d effects beyond 
the Gross- Pijtaevskii approximationtirllj and reduced di- 
mensionality!^. r~in 

As has been shown previouslyt2Hlil, periodic station- 
ary states on a ring can be characterized as bright or 
dark soliton trains for attractive or repulsive nonlinear- 
ity, respectively. It follows that our numerical studies can 
in part be understood by combining the extensive ana- 
lytical results from the physical and mathematical liter- 
ature on the stability of single solitondiHO and the na- 



ture of two-soliton interactionst^~E3 and extrapolating to 
soliton- trains. However, because we work with bounded 
systems it is necessary to consider scaling issues. Adding 
an external potential term to the NLSE adds an addi- 
tional parameter which cannot be scaled into space, time, 
or normalization of the amplitude. This extra parameter 
can be incorporated as the coefficient of the nonlinearity. 
In the context of the BEC, it has a simple physical inter- 
pretation as being proportional to the number of atoms 
in the condensate. Equivalent ly it can be incorporated 
as the coefficient of the dispersion, in which case it has 
a physical interpretation as £ 2 , where £ is the healing 
length which measures tlmdength scale of changes in the 
condensate wavefunctionEl As a result there are three 
different scale regimes: overlapping, adjacent, and well- 
separated. These terms, whic h ap ply to the spacing of 
solitons, will be defined in Sec. II A. The stability proper- 



ties and perturbation-induced dynamics of soliton-trains 
in these three regimes can differ greatly. 

The manner in which periodic soliton-train stationary- 
states on a ring respond to stochastic white noise and 
higher order nonlinearity carries directly over to the har- 
monic potential. We present a representative set of sta- 
tionary states for the harmonic potential and show how 
they can be manipulated via simple phase profiles, or 
phase engineered. It is also shown how the form of the 
ground state depends on the strength of the nonlinear- 
ity, i.e. the number of atoms. As single dark or grey 
solitons have recently been obserpjcUfor the first time in 
the BEC in a harmonic potentialc3t3, our investigations 
are of immediate benefit to BEC experiments. 

We especially focus on the physical meaning of differ- 
ent types of stability. As the NLSE has no dissipation, 
solutions to it are at best stable, but not asymptotically 
sco. Furthermore, stability on the time scale of an ex- 
periment is not the same as mathematical stability. For 
repulsive nonlinearity stationary states are generally sta- 
ble; for attractive nonlinearity, th ey axe quasi-periodically 
stable, as shall be defined in Sec. [HB 2 . In the right pa- 
rameter regimes such states can be both observable and 
useful for experimentalists. 



II. PERIODIC SOLUTIONS ON THE RING 

The analytic description of all periodic stationary 
states of the ID NLSE on a ring of circumference L is 
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reviewed. The dynamics in response to perturbation by 
stochastic white noise and higher order nonlinearity are 
then illustrated and discussed. 



A. Review of analytic solutions 

The NLSE in ID has a number of special properties 
which arepdescribed in the mathematical literature. It is 
integrableEjEj, imuz-_be solved exactly by the inverse scat- 
tering transformt3'E2l andiias a countably infinite number 
of conserved quantitiesEZlEi We have taken advantage of 
its special properties in finding its stationary solutions. 
It models the BEC in the quasi-one-dimensional regime, 
which is obtained when the transverse dimensions of the 
condensate are on the order of its healing length and its 
longitudinal dimension is much longer than its transverse 
onesE3. 

We begin with the cubic NLSE with an additional scal- 
ing factor which incorporates the strength of the nonlin- 
earity: 



[-dl ±rH/0M)l 2 ] 
[-^ 2 ±|/(x,i)| 2 ; 



f{x,t) ■■ 
f(x,t) 



id t f(x, t) , 
id t f(x,t). 



(la) 
(lb) 



The two forms of Eq. 1 have different physical emphases. 
In the first, r\ is proportional to the number of condensed 
atoms N and the s-wave scattering length between atoms; 
in the second, £ is the healing lengthO, which gives the 
length scale of variations in the condensate wavefunction 
f{x,t). As r] oc £ -2 , the two forms are equivalent. The 
± refers to repulsive or attractive two-body atomic inter- 
actions, respectively, and x and t are rescaled space and 
time. 

Assuming the wavefunction to be of the form f(x, t) = 
f(x)e~ lfJ ' t results in the stationary NLSE 



-d 2 x ± V \f(x)\ 2 ]f(x)= t xf(x), 



(2) 




where ii is the eigenvalue which corresponds to the chem- 
Dotential of the BEC. All stationary solutions to 
may be written in terms of Jacobian elliptic func- 
. B)A (properties of such functions are reviewed 
elsewherelljO'L^I. There are five normalizable, symmetry- 
breaking, periodic solution-types to the stationary NLSE 
on the ring. They are pictured in Fig. 0. All five solution- 
types are found by solving Eq. (||) subject to normaliza- 
tion and boundary conditions, and are described in detail 
in references [10] and [11]. 

Real stationary states in one-to-one correspondence 
with those of the particle-on-a-ring problem in linear 
quantum mechanics are 



f(x) = Asn 



2jK(m)x 
L 



rn 



and 



m)x 



(3) 



(4) 



for repulsive and attractive nonlinearity, respectively. 
They shall be designated sn-type and cn-type, in accord 
with their functional form. For the case where j = 2, 
Eq. (|) is shown in Fig. [jja) and Eq. (|) in Fig. 0(c). A 
is the amplitude; j G {2, 4, 6, • • •} is the number of nodes; 
K(m) is the complete Jacobian elliptic integral which 
is the quarter period of these functions; S is an arbitrary 
translational offset, which leads_to a Kosterlitz-Thouless- 
type entropy in one dimensioned; L is the circumference 
of the ring; and < m < 1 is the Jacobian elliptic pa- 
rameter. The general notatipa sn(u | m) is standard for 
Jacobian elliptic functionsESE]. 

The Jacobian elliptic parameter m governs the 
strength of the nonlinearity r]. As m — > + , r] <C 1, 
sn — > sin, and cn — > cos, respectively. This is the linear, 
sinusoidal limit which reproduces the particle-on-a-ring 
stationary solutions from linear quantum mechanics. As 
m — > 1~, r] ^> 1, sn^ tanh, and cn — > sech. These are 
the dark and bright soliton, stationary— pplutions to the 
NLSE on the infinite line, respect ivelyll3'E3. This shows 
the connection between these stationary solutions and 
solitons. On the ring the j th stationary state is a 2j 
soliton-train. 

In addition to the above two solution-types there are 
three nodeless, symmetry-breaking solution-types which 
have no analogue in the particle-on-a-ring problem in lin- 
ear quantum mechanics. The first, a solution for attrac- 
tive nonlinearity, is real: 



„ t \ a i f 2jK(m)x 
f(x) = Adn ( J 2 + * 



rn 



(5) 



where j G {1,2,3, •••}. It shall be designated dn-type. 
An example is shown in Fig. |l|(c). The other two types 
are intrinsically complex. f(x) = r(x) exp(z^(x)) and for 
repulsive nonlinearity 



r(x) z 



1 - 7 dn 



/ 2jK(m)x 
{ L 



(6) 



while for attractive nonlinearity 



r(x) 2 



f2jK{m)3 



m 



7(1 -m) , 



(7) 



where in both cases the phase must be found by numer- 
ical integration from the equation 



b\x) = 



a 



r(x) 2 



(8) 



The phase and amplitude for j = 1 and j = 3 are shown 
in Fig. 0(b) and 0(d). In the repulsive case A 2 j is the 
depth of the density minima below the constant back- 
ground. When 7 = 1 Eq. (||) is recovered. In the attrac- 
tive case 7 interpolates between the cn-type and dn-type 
solutions in Eqs. (0) and (||). In both cases < 7 < 1; 
a is a constant of integration; j G {1,2,3, •••} is the 
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number of density minima or maxima, respectively; and 
each stationary state has a complex-conjugate, degener- 
ate partner. 

The intrinsically complex solutions for repulsive non- 
linearity are interpreted as density-notch solitons mov- 
ing with speed c on the ring with an opposing momen- 
tum boost of the condensate of speed — c, which results 
in a stationary state in the lab frame. Density-notch 
solitons have a speed between zero and the Boguliubov 
sound speed, ranging from maximal to zero depth, re- 
spectivelyLd Those not of maximal depth are called grey 
solitons, while those which are of maximal depth and 
therefore form a node are called dark solitons. Fig. |l|(b) 
shows the bounded, quantized version of a stationary 
grey soliton on a ring. 

All attractive symmetry-breaking, longitudinally pe- 
riodic, stationary solution-types, i.e. the cn and dn- 
types shown in Fig. E|(b) and the intrinsically complex 
one shown in Fig. [j](d), are described by the Cj point 
symmetry group, where j is the number of peaks. There 
are j nearly degenerate solutions. For even j there is a 
real dn-type-cn-type pair and (j — 2)/2 degenerate intrin- 
sically complex pairs. For odd j there is a real dn-type 
solution and (j — l)/2 degenerate intrinsically complex 
pairs. Thus by group theory these three solution-types 
form the complete set of periodic stationary states made 
of evenly spaced peaks. 

existence of nodeless solutions requires a minimum 
Based on this bound, the effective spatial extent 
of a soliton is 2TT^/§^ or 27rv // 2£ for repulsive or attrac- 
tive nonlinearity. In our plots we shall choose £ = 1 
and let the length of the ring determine the strength of 
the nonlinearity, since on a finite interval the important 
quantity is £/L. The three scale regimes mentioned in 
the introduction are then defined as follows: j£/L <C 1 
is the well-separated regime; j£/L ~ 1 is the adjacent 
regime; and j^/L^> 1 is the overlapping regime. 

B. Dynamics under perturbation 

Of experimental interest is the stability of the above- 
mentioned stationary states on physical time scales. This 
is investigated numerically in this and the following 
sections. The algorithm employed uses a fourth or- 
der Runge-Kutta in time and a filtered pseudospectral 
method in space to propagate solutions over space-time 
intervals of interest. Although some solutions are illus- 
trated over short time scales for presentation purposes, 
all simulations were performed over time scales longer 
than experimental lifetimes of the BECEJ. In this sec- 
tion white noise and higher order nonlinearity are con- 
sidered; phase engineering of periodic solutions on the 
ring has been investigated elsewheretl. Throughout our 
presentation the level of initial white noise will represent 
a number of uncondensed atoms on the order of 1% of 
the condensed onesS 



Kolomeisky et alxB showed that for the BEC a true 
one-dimensional mean field theory, rather than the quasi- 
one-dimensional regime of the three-dimensional NLSE, 
has quintic rather than cubic nonlinearity. For this rea- 
son a quintic perturbation is chosen for higher order non- 
linearity. Transverse confinement on the order of the 
healing length should be described by the NLSE, while 
transverse confinement on the order of the atomic inter- 
action length should be described by a one-dimensional 
mean field theory. As these scales are two or three orders 
of magnitude apart, the combination of cubic and quintic 
nonlinearity may describe the intermediate regime. This 
full range of confinement scales-. should soon be available 
for experimental investigationEA 

1. Repulsive nonlinearity 

In response to white noise, density-notch soliton-trains 
drift but are otherwise stable. This may be understood 
by extrapolation from the stability properties of single 
solitons and two-soliton interactions. Density notch soli- 
tons repel elastically. Single solitons respond to stochas- 
tic perturbation by emission of radiation and by a change 
in velocity and depth. Velocity and depth are a function 
of a single parameter, sometimes called the soliton phase 
angle. A full analysis requires separate consideration of 
the constant background and the density notch. Kivshar 
and YangfJ have completed this analysis for single soli- 
tons on the infinite line by using a variational approach. 

In Fig. ||(a)-(b) the evolution of a sn-type station- 
ary state described by Eq. (||) with initial white noise 
is shown. This is the lowest such state on the ring. On a 
longer time scale there is some drift. In this case the two 
density notches retain their spacing and phase relation- 
ship, which shows that, as with single solitons, only the 
velocity of the soliton train changes. In the box this same 
stationary state is the first excitation above the ground 
state. Numerical studies, though not shown here, reveal 
the single density notch in the center drifts randomly 
between the box boundaries until it comes close to one, 
at which point it reverses direction abruptly. This can 
be explained by interpreting the falling off of the wave- 
function at the boundaries as pinned solitonsL3, which 
repel the central density notch when it comes within the 
interaction length of 7r\/6£. 

In Fig. ||(c)-(d) the evolution of an intrinsically com- 
plex stationary state described by Eqs. (||) and (||) with 
initial white noise is shown. This is the lowest symmetry- 
breaking state on the ring. It consists of both a back- 
ground phase ramp and a density notch, and is quite 
stable. Because phase quantum number is conserved for 
the NLSE, the background cannot be perturbed into a 
lower energy state. And a grey density notch has the 
same response to perturbation as a dark density notch. 
Therefore the two together emit radiation and, on longer 
time scales, drift, but remain otherwise unchanged. 
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In Fig. ||(a) the evolution of a highly excited state is 
shown. This demonstrates that for repulsive nonlinear- 
ity there is no difference between the stability proper- 
ties of soliton-train stationary states in the overlapping, 
adjacent, or well-separated regimes, a fact not a priori 
apparent from single soliton considerations. In Fig. ||(b) 
continuous white noise is added. This perturbs the equa- 
tion rather than the initial solution. As there is no dissi- 
pative term in the NLSE, the noise builds up over 1000 
time units to a very high level. Even in this case, the 
only instability exhibited by a highly excited, intrinsi- 
cally complex solution is drift. 

Finally, in Fig. ^(c)-(d) quintic nonlinearity is added 
to the same sn-type stationary state as was used in 
Fig. ||(a)-(b). Eq. (1) then becomes: 

[-dl ± r, | f{x, t)\ 2 +b\ f{x, t) | 4 ] f{x) = tdtfix, t) (9) 

We used b = ±0.1 for (c) and (d), respectively. In both 
cases the initial state was found to be stable with respect 
to this perturbation, with only small modulations in the 
background and a slight reshaping of the density notchH. 
For intrinsically complex stationary states a similiar set 
of results was obtained, except that the balance between 
the constant background and the grey density notch was 
no longer the same, so that the solution translated. This 
suggests that soliton-train stationary states may exist for 
transverse confinement of the BEC on scales intermediate 
between the atomic interaction length and the healing 
length. 

2. Attractive nonlinearity 

The study of soliton trains differs from that of single 
solitons or two-soliton interactions because it no longer 
suffices to consider the well-separated limit. The exact 
equation for interactions betwepn two bright solitons has 
been derived in full by GordonE3. As it is difficult to in- 
terpret when the solitons are not resolvable, he has con- 
sidered the non-overlapping limit and found 

q = -4exp(-2g) cos(2^) 
^ = 4 exp(-2<?) sin(2^) (10) 

where 2q is the separation of the two solitons and 2\I/ is 
their relative phase. The interaction depends exponen- 
tially on the separation-and sinusoidally on the relative 
phase. Desem and ChuLJ used this work to evaluate in- 
teraction minimization schemes in the context of fiber 
optics. 

We found that the key to understanding the stabil- 
ity of intrinsically complex stationary states perturbed 
by initial white noise for attractive nonlinearity lay in 
the overlapping regime of these two soliton interactions. 
Figure [| shows the results of our studies. Following the 
six subpanels of the figure from left to right and top to 



bottom, it is apparent that the topology of the space- 
time profile changes from connected to disconnected as 
the phase varies from to n. A hole opens in the cen- 
ter, rotates, and then separates the density profiles of the 
two solitons. Thus the interaction changes continuously 
from completely attractive to completely repulsive. For 
2\I> G [7r, 2tt] the topology reverses itself in the same way. 

Single bright solitons respond to stochastic perturba- 
tion not only by emission of radiation and by a change 
in velocity, but also by a drift in phaseO. As a soliton 
train requires a fixed phase relationship between its com- 
ponents, it is to be expected that the stability properties 
of solitons trains for attractive nonlinearity differ from 
those of repulsive nonlinearity, as may be seen for exam- 
ple hythe numerical studies of densely-packed solitons by 
ArbeM 

In Fig. ||(a)-(b) the evolution of cn-type stationary 
state described by Eq. (||) is shown with initial white 
noise. In the adjacent regime this solution is quite stable 
out to very long time scales. However, in the overlapping 
regime such a stationary state undergoes interactions, 
while in the well-separated regime the peaks behave as 
individual solitons, as may be seen in Fig. |(c)-(d). The 
drift in phase is especially evident in (d), where in the 
first interaction there is clearly density exhange between 
the two solitons, despite their initial phase difference be- 
ing 7T. 

In Fig. ||(a)-(b) the evolution of a real dn-type sta- 
tionary state described by Eq. (^|) is shown with initial 
white noise in the adjacent regime. As the phase differ- 
ence between peaks is zero, it is not surprising that the 
solution seems to quickly go unstable. However, there 
are quasi-periodic recurrences, which are noted with as- 
terisks in (b). For this reason the solutions are termed 
quasi-periodically stable. This is especially evident on 
the longer time scale shown in Fig. ||(c), where it may be 
seen that though the peaks continue to exchange mass 
the overall integrity of the soliton-train remains intact. 
The drift of the train is similiar to that found for the 
cn-type stationary states in the adjacent regime, and as 
before, in the overlapping regime it is unstable while in 
the well-separated regime the peaks behave as indepen- 
dent solitons. 

The intrinsically complex stationary states have sim- 
iliar properties to the real dn-type stationary states: in 
the adjacent regime they evolve and retain their overall 
integrity, and in other regimes they are unstable, as is 
shown in Fig. @(a)- (b). However, even in the unstable 
case they do not ever superimpose to make a sharp peak. 
This is due to the conservation of phase quantum number 
by the NLSE. For the BEC, this is an important point. In 
higher dimensionality solutions to the NLSE with attrac- 
tive nonlinearity collapse. This occurs when the density 
becomes large enough for the nonlinear term to domi- 
nate the kinetic energy. A highly excited, intrinsically 
complex stationary state, though unstable, can lower the 
maximum density by a factor of forty or more, thereby 
preventing collapse. 
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Finally quintic nonlinearity is added, as in Eq. (||). 
For b = 0.1 the cn-type stationary state in the adjacent 
regime pulses but is otherwise unaffected, as may be seen 
in Fig. 0(c). In Fig. 0(d) it is seen that for b = —0.1 
the same stationary state undergoes large oscillations but 
does not lose its integrity as a soliton-train. Although not 
illustrated here, the ground state at this scale, a single 
density peak, deforms slightly but is otherwise unaffected 
for b = ±0.1. 



III. HARMONIC POTENTIAL 

Since many BEC experiments involve magneto-optical 
trapsO, the consideration of the stability properties of 
the stationary states of such a potential has immediate 
experimental relevance. A harmonic confining potential 
is added to Eq. (1) and it is rescaled to the form 

[-d 2 x ± T) | /(*, t) | 2 W ] /(*, t) = idtffa f), (11) 

where a measures the strength of the harmonic trap and rj 
the strength of the nonlinearity, which is proportional to 
the number of condensed atoms. As before, we consider 
the class of stationary solutions to Eq. (JTTj) by letting 

f(x,t) = f(x)e~ %flt , so that 

\-dl ± t? | f(x, t) | 2 +ax 2 ] f[x) = (12) 
is the resulting eigenvalue problem with eigenvalue \i. 



A. Stability of stationary states 

The addition of the harmonic potential no longer al- 
lows closed form solutions to be obtained. However, the 
normalized eigenfunctions of Eq. (|l^) can be constructed 
numerically via standard shooting methods&3. Fig. || 
depicts the stationary solutions of Eq. (|l^ ) normalized 
to unity with a = 0.02 and r\ = 1. In Fig. |](a)-(b) 
the zeroth, first, and sixth modes are depicted for at- 
tractive and repulsive nonlinearity, respectively. In the 
lower modes it is seen that attractive nonlinearity sharp- 
ens the peaks and troughs while repulsive nonlinearity 
makes them spread out. Highly excited states, as for ex- 
ample the sixth mode, are further in the linear regime, 
and therefore resemble the Hermite polynomial station- 
ary solutions to the analogous problem in linear quantum 
mechanics. These results are qualitatively in accord with 
the sn-type and cn-type closed formtj'EiJ stationary states 
for box or periodic boundary conditions. 

As the condensate is grown the effective number of 
atoms is increased, thereby increasing the strength of 
the nonlinearity via the parameter 77. To illustrate this, 
Fig. ^| depicts the density of the zeroth mode as a func- 
tion of increasing nonlinearity. Here 77 varies from to 
3 and the eigenfunctions are clearly seen to deform as 77 
increases. For attractive nonlinearity, the ground state 



steepens significantly from the original Gauss-Hermite 
polynomial which is the solution for the linear case of 
77 = 0. In contrast, repulsive nonlinearity causes the 
Gauss-Hermite polynomial to spread out as the nonlin- 
earity is increased. It takes on a parabolic shape with 
Gaussian tails, known as the Thomas-Fermi limit for the 
BEC. Thus as atoms are condensed in a harmonic poten- 
tial trap, the effect of the mean-field contribution is to 
significantly reshape the stationary states. 

To study the stability of these states a substantial 
initial white-noise perturbation was added and then 
Eq. dll| ) was solved. The zeroth through sixth modes 
were tested, but for the sake of brevity only the first and 
sixth ones are illustrated. In accord with the results of 
Sec. [n] it was expected that all solutions would be stable 
for repulsive nonlinearity, while for attractive nonlinear- 
ity they would be stable in the adjacent but not the over- 
lapping regime. Note that here there is no well-separated 
regime, due to the effect of the harmonic potential. 

Fig.[H](a)-(b) and (c)-(d) depicts the results of the sim- 
ulations for a = 0.02 and 77 = 1 for a time of t = 1000 for 
attractive and repulsive nonlinearity, respectively. Both 
the perturbed attractive and repulsive stationary states 
are observed to be stable over very long times. Fig- 
ure |l(](b) however, starts to exhibit an oscillatory in- 
stability at long times, much like that observed for the 



overlapping regime in Sec. [IB 2 



Although not illustrated here, the effect of either nega- 
tive or positive quintic n onlinea rity of the form and mag- 
nitude described in Sec. II B 1 is to slightly deform the 
eigenstates but leave them otherwise unaffected. 



B. Phase engineering 

Phase engineering has been used to successfully cpate 
solitons in BECs confined in a harmonic potentiaHLj. 
Specifically, a step function in the phase is used to create 
a density notch which then propagates across the conden- 
sate. It has been shown elsewhereoH that in quasi-one- 
dimensional confinement both bright and dark solitons 
may be manipulated, or phase engineered, by various 
simple phase profiles. Here it is shown that the same 
is true of excited states in a harmonic potential. 

To induce dynamics in the stationary states, the eigen- 
function solutions depicted in Fig. || were modified by in- 
troducing the following two phase profiles into the initial 
conditions: 



I. f(x,0) = f(x)exp(ipx) 

II. f(x,0)=f(x)exp(il3\x\), 



(13) 
(14) 



where /3 determines the phase pertubation strength. For 
profile I, all peaks are ramped in the same fashion, 
whereas for profile II, the peaks are ramped in opposite 
directions initially depending upon their location in x. 

Fig. |ll] illustrates the resulting dynamics for both at- 
tractive and repulsive nonlinearity. In (a) and (c), phase 
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profile II is imposed on the first mode for a = 0.02, r\ = 1, 
and (3 — 0.3. This initially leads to a repulsion of the 
peaks. However, the potential counteracts this effect and 
the peaks undergo an oscillatory particle-like motion. In 
contrast, phase profile I keeps the peaks moving in unison 
and the potential once again acts to trap the peaks in- 
side the potential by generating a periodic motion within 
the well. These two cases are analogous to the oscillatory 
eigenmodes of the coupled pendulum, and clearly demon- 
strates the soliton-like behaviour of the stationary states 
when phase engineered. Thus in the mean- field approx- 
imation the BEC, itself consisting of particles, has solu- 
tions of a particle-like nature. 



IV. CONCLUSION 

The stability of excited stationary states of a 
trapped gaseous Bose-Einstein condensate in a quasi-one- 
dimensional environment has been studied with respect 
to addition of stochastic white noise, externally imposed 
phase profiles (phase engineering), and perturbations of 
the effective Hamiltonian by inclusion of interactions pro- 
portional to higher powers of the density. This study has 
been carried out within the context of investigating per- 
turbations to the simplest mean-field picture: that which 
gives rise to the cubic nonlinear Schrodinger equation. 
Such excited states are soliton-like, and, as such, 
recently been observed experimentally for the BEClii 

Two topologies were considered: periodic solutions on 
a ring and confinement in a harmonic potential. It was 
found that in both physical contexts all of the above 
types of perturbations basically left the solitonic struc- 
tures intact, even when the initial stationary states had 
become unstable. In quasi-one-dimension, harmonic con- 
finement enhanced stability. For attractive nonlinear- 
ity the perturbed stationary states were either stable or 
quasi-periodically stable. The most stable solution-type 
was one for which alternating peaks were-placed adjacent 
to one another and separated by nodea23'E3. For repul- 
sive nonlinearity, excited states drifted but were other- 
wise stable. 
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FIG. 1. Shown are the five periodic stationary solu- 
tion-types in one dimension, (a) Amplitude of real solution 
and (b) amplitude and phase of intrinsically complex solu- 
tion for repulsive nonlinearity; (c) amplitude of dn-type and 
cn-type solutions and (d) amplitude and phase of intrinsically 
complex solution for attractive nonlinearity. 
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FIG. 2. Stability for repulsive nonlinearity, periodic solu- 
tions on the ring, (a) Density and (b) phase of real, sn-type 
solution with initial white noise; (c) density and (d) phase of 
intrinsically complex solution. The lowest energy solution is 
shown for each type. A short time scale was used to illumi- 
nate the phase but the same stability properties hold over an 
order of magnitude longer time. The phase is plotted mod 
2tt. 
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FIG. 3. More stability properties for repulsive nonlinear- 
ity. (a) Highly excited intrinsically complex solutions are sta- 
ble to perturbation by initial white noise, despite slow drift, 
(b) Continuous noise builds up because there is no dissipation 
term in the NLSE, but even highly excited solutions remain 
stable. For (c) positive and (d) negative quintic nonlinear- 
ity without white noise the stationary states develop slight 
fluctuations but remain stable on very long time scales. 
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FIG. 4. Phase dependence of bright soliton interactions. 
Shown are the space-time projection of the density of colliding 
bright solitons. The phase difference between the two solitons 
is (a) 0, (b) tt/16, (c) tt/8, (d) tt/4, (e) tt/2, and (f) tt. The 
time scale in each case is 50 time units and the interaction 
length is 2 r K\[ 7 i. Note that the circumference of the ring is 
much greater than the length scale shown here. 




FIG. 5. Stability and the importance of scale, attractive 
nonlinearity. (a) Density and (b) phase of the cn-type so- 
lution with initial white noise in the adjacent regime. The 
phase is plotted mod 2tv. However, in the (c) overlapping and 
(d) well- separated regimes the perturbed cn-type solution is 
unstable. 




FIG. 7. More stability properties for attractive nonlin- 
earity. (a) The long time behaviour of intrinsically complex 
solutions in the adjacent regime when perturbed by initial 
white noise is the same as for dn-type solutions, (b) In the 
overlapping regime, intrinsically complex solutions are unsta- 
ble but do not collapse to a ground state, due to conserva- 
tion of phase quantum number. Positive quint ic nonlinearity 
causes the peaks to pulse while (d) negative quintic nonlin- 
earity leads to large oscillations but not actual breakup of the 
solitons. 




FIG. 6. Quasi-periodic stability and the importance of 
scale, attractive nonlinearity. (a) Density and (b) phase of 
the nodeless, dn-type solution with initial white noise in the 
adjacent regime. The asterisks in (b) mark the recurrence of 
the solution, clearly visible here on a short time scale, (c) 
On longer time scales but still in the adjacent regime, the 
peaks exchange mass and drift, but continue to recur, (d) In 
the well-separated regime the density peaks behave as inde- 
pendent solitons. Although not shown here, the perturbed 
dn-type solution is also unstable in the overlapping regime. 
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FIG. 8. Stationary states of the harmonic potential for 
(a) attractive and (b) repulsive nonlinearity. Pictured are the 
ground state, first mode, and sixth mode. Note that states in 
the overlapping regime, such as the sixth mode, are predomi- 
nantly linear, while those in the adjacent regime, such as the 
first mode, are strongly nonlinear. There is no well-separated 
regime in the harmonic potential. 
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FIG. 9. Dependence of the ground state on the strength of 
the nonlinearity for the harmonic potential, (a) For attractive 
nonlinearity the ground state deviates rapidly from the linear 
gaussian solution towards a sharp peak, (b) For repulsive 
nonlinearity it broadens to become parabolic with gaussian 
tails. 
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FIG. 10. Stability of stationary states in the harmonic 
potential. The first and sixth modes are propagated with 
initial white noise over long time scales for (a)-(b) attractive 
nonlinearity and (c)-(d) repulsive nonlinearity. The stability 
properties are identical to those of the analogous solutions for 
periodic solutions on the ring: only (b) shows instability, as 
it is for attractive nonlinearity in the overlapping regime. 




FIG. 11. Phase engineering of stationary states in the har- 
monic potential. The particle-like nature of these solutions is 
exhibited for both (a)-(b) attractive nonlinearity and (c)-(d) 
repulsive nonlinearity. In (a) and (c) an equal and opposite 
linear phase ramp was used on the two density peaks, while 
in (b) and (d) an identical phase ramp was used. This pro- 
duces two oscillatory modes analogous to those of the coupled 
pendulum. 
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